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ABSTRACT 

A spherical Lagrangian hydrodynamical code has been written to study the formation 
of cosmological structures in the early Universe. In this code we take into account 
the presence of collisionless non-baryonic cold dark matter (CDM), the cosmological 
constant and a series of physical processes present during and after the recombination 
era, such as photon drag resulting from the cosmic background radiation and hydrogen 
molecular production. We follow the evolution of the structure since the recombination 
era until the present epoch. As an application of this code we study the formation of 
voids starting from negative density perturbations which evolved during and after 
the recombination era. We analyse a set of COiJi^normalized models, using different 
spectra to see their influence on the formation of voids. Our results show that large 
voids with diameters ranging from 10/i _1 Mpc up to 50ft, _1 Mpc can be formed in a 
universe model dominated by the cosmological constant (J7a ~ 0.8). This particular 
scenario is capable of forming large and deep empty regions (with density contrasts 
5 < —0.6). Our results also show that the physical processes acting on the baryonic 
matter produce a transition region where the radius of the dark matter component is 
greater than the baryonic void radius. The thickness of this transition region ranges 
from about tens of kiloparsecs up to a few megaparsecs, depending on the spectrum 
considered. Putative objects formed near voids and within the transition region would 
have a different amount of baryonic/dark matter when compared with Ob/fid- If one 
were to use these galaxies to determine, by dynamical effects or other techniques, the 
quantity of dark matter present in the Universe, the result obtained would be only 
local and not representative of the Universe as a whole. 

Key words: hydrodynamics - galaxies: formation - cosmology: observations - cos- 
mology: theory - dark matter - large-scale structure of Universe. 



1 INTRODUCTION 

In the 1970s a series of studies showed the existence of fil- 
amentary distribution of galaxies in clusters and superclus- 
ters and large void regions among such filaments. Gregory 
& Thompson (1978) found in galaxy redshift surveys, in the 
direction of the Coma cluster, a filamentary structure that 
is part of the nowadays denominated 'Great Wall' (Geller & 
Huchra 1989). This supercluster is encompassed by a large 
region of low mass density (or void). 

Additional investigations by Gregory, Thompsom & 
Tifft (1981) and Chincarini, Rood & Thompson (1981) 
showed similar results for the Hercules and Perseus-Pisces 
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superclusters and, later, other redshift surveys showed that 
the superclusters are boundaries of underdense regions (Gio- 
vanelli, Haynes & Chincarini 1986; de Lapparent, Geller & 
Huchra 1986; da Costa et al. 1988). 

A redshift survey in the Bootes region (Kirshner et al. 
1981) discovered a large void structure whose estimated di- 
ameter is ~ 60/i~ 1 Mpc. Later, observations showed that this 
region in Bootes is not totally empty of matter - it contains 
some galaxies (Thuan, Gott & Schneider 1987; Dey, Strauss 
& Huchra 1990). In spite of this, the void in Bootes has a 
lower mass density as compared to the mass density of the 
Universe. 

In particular, Dey et al. (1990) found 21 galaxies (of 
which 13 are IRAS sources) in Bootes, which led the au- 
thors to estimate the average density contrast of this re- 
gion to be —0.84 < 5 < —0.66. More recently, Szomoru et 
al. (1996a,b) analysed, through the VLA (Very Large Ar- 
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ray) technique, ~ 1 per cent of the volume of the void in 
Bootes (~ 1100/t _3 Mpc 3 ). These authors found that galax- 
ies in Bootes are systems rich in gas, composed by late-type 
galaxies. The study of their optical properties shows that 
these systems are very similar to field galaxies of the same 
morphological type. 

From the theoretical point of view, some models of 
galaxy formation (Dekel & Silk 1986) that studied the cold 
dark matter (CDM) scenarios including a bias parameter, 
predict that dwarf galaxies should originate from the la 
fluctuations, thus being more smoothly distributed than the 
rare high-density peaks that form the more massive galax- 
ies. In these scenarios, the dwarf galaxies could trace the 
dark matter and also fill up the voids. On the other hand, 
Popescu, Hopp & Elsasser (1997), using results of a search 
for emission-line galaxies (ELGs) , did not find that the voids 
are occupied by a homogeneous population of dwarf galax- 
ies. They found in their sample a few galaxies near the void 
regions, but the number of void galaxies is not significant as 
compared to the field galaxies. 

It is worth stressing that, although the Bootes void has 
been one of the first underdense region studied, the distri- 
bution of galaxy clusters present in the Abell catalogue also 
reveals the existence of large empty regions with diameters 
of ~ (20 - 50)h -1 Mpc (Bahcall & Soneira 1983; Huchra et 
al. 1990). 

Some authors (Goldwirth, da Costa & Van de Wey- 
gaert 1995) argue that the typical size of underdense regions 
(voids) is related to the first zero of the correlation function 
£(r), which gives ~ 20/i _1 Mpc. This length is less than the 
void region sizes found in studies of the large-scale structure 
of the Universe. 

Recently, El- Ad, Piran & da Costa (1996), using a par- 
ticular way to define voids, have found an average diameter 
of 37 ± 8/i _1 Mpc. In their particular method to define voids 
they allow the existence of field galaxies therein, this be- 
ing a very reasonable assumption because galaxies could in 
principle be formed therein or escape to there, as a result of 
mutual galaxy interactions. 

Also, quite recently, Miiller et al. (2000) studied the dis- 
tribution of void regions in the two-dimensional slices of the 
Las Campanas Redshift Survey (LCRS; the deepest redshift 
survey presently available). They found that the distribution 
of void sizes scales with the mean galaxy separation, A. In 
particular, they found that the size of voids covering half of 
the area is given by D mcd ~ A + (12±3) = (20-30)h _1 Mpc, 
where D me d is the median void diameter. 

In the past, the models addressed to explain the forma- 
tion of voids in the distribution of galaxies considered their 
formation as a result of the process of clustering (Aarseth 
& Saslaw 1982; Peebles 1982; Vettolani et al. 1985) or as 
a result of the evolution, after the recombination era, of 
negative density perturbations (Hausman, Olson & Roth 
1983; Hoffman & Shaham 1983; Filmore & Goldreich 1984; 
Bertschinger 1985; de Araujo & Opher 1990, 1993). 

The formation of voids as a result of explosions of pri- 
mordial objects, such as quasars or other pre-galactic objets 
(Ikeuchi, Tomisaka & Ostriker 1983) was also considered. 
This scenario, however, might have some problems because 
of the distortions that could be caused on the cosmic back- 
ground radiation. On the other hand, this explosive scenario 
was revisited by Miranda & Opher (1996, 1997) and the for- 



mation of voids through such a mechanism may be consis- 
tent, within some limits, with the C'OBE satellite data. 

In an interesting study, Srianand (1997) has reported 
the detection of a large void ~ 7h~ Mpc centred near the 
quasi-stellar object (QSO) Tol 1037-2704, where this void 
could be produced by an excess of ionization due to the 
QSO, thereby constituting another mechanism to produce 
voids of several megaparsecs. 

Also, de Araujo & Opher (1997) considered an unusual 
possibility for the formation of void regions, namely, through 
putative primordial magnetic fields. 

However, after many years of studies it is not completely 
clear how voids are really formed, particularly the largest 
ones, such as those found by El- Ad et al. (1996) and Miiller 
et al. (2000). It is possible to have their formation through 
any of the mechanisms mentioned above. 

It is worth mentioning that, although the dimensions 
of voids depend on the particular technique employed to 
find them (see e.g. Lindner et al. 1995, 1999; El-Ad et al. 
1996; EL-Ad & Piran 1997; Miiller et al. 2000), their very 
existence is widely accepted. 

Concerning the presence of non-baryonic dark matter 
in the structures of the Universe, there is considerable evi- 
dence for its existence. It has been argued that non-baryonic 
dark matter can dominate the density and the dynamics of 
the Universe (for reviews see Kolb & Turner 1990; Primack 
1997). Thus, there is no reason to believe that dark mat- 
ter is not present inside voids; therefore, it is important to 
take into account this kind of matter in the study of the 
formation of underdense regions, in order to verify its role. 

It is also important to stress that a calculation taking 
into account the non-baryonic dark matter and a series of 
effects during and after the recombination era, focused on 
void formation, has so far not been performed. 

In the present work we use a spherical Lagrangian hy- 
drodynamical code, written by one of us (ODM), which al- 
lowed us to study the formation of cosmological structures. 
Besides the presence of the non-baryonic dark matter, we 
consider in our calculations the expansion of the Universe, 
the photon drag due to the cosmic background radiation, 
the recombination processes and molecular hydrogen forma- 
tion. We also included, in our model, the presence of the 
cosmological constant to see its influence on the dimension 
of void regions. 

In Section 2 we describe the model, in Section 3 we 
present the numerical results and the discussions, and in 
Section 4 we present the conclusions. 



2 THE MODEL 

We assume spherically symmetric negative density pertur- 
bations which produce clouds of baryonic and non-baryonic 
dark matter with densities lower than the density of the 
Universe. We treat the dark matter and the baryonic mat- 
ter as two fluids coupled by gravity. We assume, as an initial 
condition, a top-hat density profile for the void region. 

The hydrodynamic equations that describe the dynam- 
ics of the baryonic density perturbations are 

-dt+-,Q- r ^ *>«b)=0, (1) 
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(the mass conservation equation to the baryonic matter writ- 
ten in spherical coordinates), where pb is the baryonic mass 
density, Vb is the velocity and r is the radial coordinate. 
The equation of motion is given by 



Dvb 
Dt 



1 dP 
Pb Or 



4 a T bT?x a 
3 m D c 



H(t)r] 



(2) 



where P is the pressure of the gas, x e is the degree of ioniza- 
tion, b = 4<tsb/c (with a sb the Stefan-Boltzmann constant), 
T r is the temperature of the radiation, m p is the proton 
mass, c is the velocity of light and H(t)r is the Hubble flow, 
which takes into account the expansion of the Universe. The 
last term on the right-hand side of equation (2) describes the 
photon drag and D/Dt = d/dt + vd/dr is the convective 
derivative. 

The field equation (Poisson's equation) is 



V 2 $ = 47rG(p b +p d )-Ac 2 



(3) 



where pd is the mass density of the non-baryonic dark com- 
ponent and A is the cosmological constant. 

The density parameters of the Universe read 



where 



Ac 2 
3H* 



and Qk 



kc 2 



(4) 



(5) 



In the equations (4) and (5), fib (fid) is the bary- 
onic (non-baryonic dark matter) density parameter, 7? is the 
scalefactor of the Universe, Ho is the Hubble parameter (the 
zero subscript stands for its present value) and k is the cur- 
vature of the Universe, that is, ft = — 1 for an open universe, 
k — for a flat universe and k = +1 for a closed universe 
(thus fik is the curvature parameter). 

The energy equation is written as 

DE = P Dpb 
Dt p\ Dt 

where E is the thermal energy per gram and L is the cooling- 
heating function of the gas. 

The cooling-heating function L is given by the summa- 
tion of four mechanisms: 



(6) 



L = Z/r + Lc + Lb_2 + L a . 

(i) The cooling from recombination Lr (see, 
Schwarz, McCray & Stein 1972) is given by: 

The ionization fraction is given by 



Dx c 

~dT 



C 



-(B 1 -B 2 )/k B T 1 



'(1 



+ 1, 



(7) 
e.g., 

(8) 
(9) 



where T m is the temperature of the matter, Na is the Avo- 
gadro number, Bi is the bound energy of the ground state 
and £>2 is the bound energy of the first excited state. The C 
and f3 parameters present in equation (9) are given by 



C = 



A 2 



A 2 



and P — a 



(27rm c fc B r,.) 3/2 
h 3 ' 



-B 2 /k B T r 



(10) 



In the above equations A2 S ,is = 8.227 s 1 is the de- 
cay rate from the 2s state to the Is state, a = 2.84 x 



10~ 11 Tm 1/ ' 2 crr^s" 1 is the recombination rate, m c is the 
electron mass and fca is the Boltzmann constant. The col- 
lisional ionization rate I (see, e.g., Defouw 1970) is given 
by 

/ = 1.23 x 10- 5 A A a;c(l + x c )p^-T^ 2 e ~ Bl/k ^ s _1 . (11) 

(ii) The Compton cooling-heating mechanism Lc ( Pee- 
bles 1968) is 



Lc 



AksNarbT^Xe 



(12) 



Lh2 = 



(iii) For the cooling by molecular hydrogen Lh2, we have 
A M 



Pb 



(13) 



The equation for Am and the complete set of equations for 
the formation and destruction of H2 molecules can be ob- 
tained in de Araujo & Opher (1988, 1989) and Oliveira et 
al. (1998a). 

(iv) The Lyman a cooling L a (Carlberg 1981) is given 

by 



L a = 1.25 x 10- n Ci2 Ml 



A27 + C21 ' ^ ^ 

where A2 7 = 8.227 s^ 1 is the two- photon emission rate, 
C21 = 1.2 x 10 _6 T,n 1 ^ 2 x nu s _1 is the collisional de- 
excitation rate, nu is the numerical density of neutral hy- 
drogen and C 12 = 2C 2 ie- 119xl ° 5 / T ». 

The equation of state, to complete the system of equa- 
tions to the baryonic matter, is 



P = k B N A Tp h (l+x c ). 



(15) 



It is worth stressing that the temperatures in our cal- 
culations are never greater than ~ 5000K; in this way the 
cooling and heating mechanisms considered here account for 
the thermal history of the baryonic matter present in the 
underdense regions that we study. 

The equations of conservation of mass and momentum 
for the dark matter are, respectively, 



dpa 1 d . 2 

— + ^ — ( r Pd«d) 



dt 

Dvd 
Dt 



dr 
= -V$. 



= 0, 



(16) 
(17) 



We consider that the negative density fluctuations 
evolve within a medium of density p t = po(l + z) 3 , where 
Po — PbuO+PduO (the subscript 'u' stands for the physical pa- 
rameters of the Universe). Thus, in our model we take into 
account the influence of the external medium on the evo- 
lution of the underdense region, and therefore as external 
boundary condition we have the pressure of the Universe, 
namely 



Pbu = fcBA r ATbuPbu(l + a;bu). 



(18) 



It is also worth stressing that in many studies of structure 
formation the influence of the external medium on the struc- 
ture under study is not taken into consideration. 

In order to calculate the temperature of the matter of 
the Universe we consider in the cooling function only the 
contributions due to the Compton cooling-heating and the 
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cooling due to the recombination processes. The other phys- 
ical processes considered for the evolution of the baryonic 
density perturbation are not relevant for the evolution of 
the matter temperature of the Universe. The energy equa- 
tion thus reads 

DE 

bu -fbu u 



Dt 



Pbu 



Dt 



— Lr — he- 



(19) 



Similarly, the collisional ionization is not important for 
the matter of the Universe. Thus, we use the study of Peebles 
(1968), where the degree of ionization is given by 



Dxbu 



/3e 



-(B 1 -B 2 )/k B T, 



Dt 
where now 
a = 2.84 x 10 



r (l - £bu) 



-1/2 



(20) 



(21) 



Equations (19) and (20) are, respectively, the equations of 
conservation of thermal energy and ionization balance, with 
Tbu the temperature of the matter of the Universe and Xbu 
its ionization fraction. 

The baryonic ambient density, i.e., the baryonic density 
of the Universe, is given by 



Pbu = Pbu()(l + z'Y 



3Hq f^b / T r 



8ttG 



where z is the redshift and To is the present cosmic back- 
ground radiation temperature. 

The equation for the time evolution of the radiation 
temperature and, therefore, the time evolution of the Hubble 
parameter is given by 



£ )=~H{t) 



-Ho 



fim(|) + (1 - On - Oa) (Jjf) +Va 



1/2 



(23) 



where f2 m = fib + fid is the density parameter of the matter. 

To solve the above set of equations, a hydrodynamical 
Lagrangian numerical code has been written, which basi- 
cally follows the method developed by Rychtmyer & Mor- 
ton (1967). It is important to stress that the total amount 
of baryonic matter and non-baryonic dark matter are taken 
to be constant throughout the calculations. 

We use in our models 9000 shells for the baryonic matter 
and 9000 shells for the dark matter. Models with 45 000 
shells for each component showed no considerable difference 
in the results. 

We also test our code, with good agreement, against 
the one-zone code of de Araujo & Opher (1991) for the col- 
lapse of baryonic and non-baryonic dark matter Population 
III objects (positive density perturbation). Also, as an ad- 
ditional test we reproduce the results of de Araujo & Opher 
(1993) for voids formed from negative density perturbations, 
where a one-zone model was used. The program conserves 
mass and total energy of the system to at least one part in 
10 6 . 



3 NUMERICAL RESULTS AND DISCUSSIONS 

In the present paper we use different power spectra of den- 
sity perturbations to study their role on the scale of voids 
formed and to see if cosmological parameters could be con- 
strained. 

First we take a typical CDM power spectrum of pertur- 
bations given by (Padmanabhan 1993) 



| 4 



Afc (4-2 Q ) 



(l + Bk + Ck 3 / 2 + Dk 2 ) 2 



(24) 



where A is the spectrum amplitude, B = 1.7(£l/i 2 ) x Mpc, 
C = 9.0(fi/i 2 )~ 3/2 Mpc 3/2 and D = 1.0(f2/i 2 )- 2 Mp C - 2 . In 
the above equation the power index is related to a by the 
relation n = (4 — 2a). The variance of fluctuations on scale 
ro reads 



AVo) = JL 



dk k 



W 2 (kr ) 



(25) 



We analysed two different window functions W(kro): 
Gaussian and top-hat filters. Our results are very weakly 
dependent on W(kro) so, for convenience and to compare 
our results with previous studies, we present here the results 
using a top-hat window function. Then, we have 



sin(fcro) cos(fcro) 



(fer ) 



(fcr ) 2 



with the volume of the top-hat window given by 



<± 3 



The mass fluctuation in a radius R reads 



— I = 2^2 / /■"'//•■|<>'i*->\ 



(26) 



(27) 



(28) 



The input parameters that describe the cosmological sce- 
narios are taken from a recent study performed by Cen 
(1998a,b), who analysed mass and autocorrelation functions 
of rich clusters of galaxies from linear density fluctuations. 
Cen presents 32 different cosmological models to determine 
the fitting parameters in his Gaussian peak method, using 
them to obtain the cluster mass function. In particular, we 
take four of these models that present different characteris- 
tics related to the present radii and present density contrasts 
of the voids. The models chosen have their parameters de- 
scribed in Table 1 (see also Cen 1998a, b). 

All models considered here have 'excess power' EP = 
1.30 ± 0.15. ('Excess power' is denned as EP = 3.4<725/as, 
where 025 is the linear rms density fluctuation in a 
25ft~ Mpc top-hat sphere at z = 0, and as has a simi- 
lar definition in an 8/i~ 1 Mpc top-hat sphere.) This param- 
eter describes the shape of the power spectrum on scales 
~ 8 — 300h~ Mpc and its values range from 1.15 to 1.45 to 
fit the COBE data. 

In the first column of Table 1 we present the number of 
the model [we use here the same numbers presented by Cen 
(1998a,b)]; in the second and third columns we present the 
density parameter for the matter and that associated with 
the cosmological constant, respectively; the fourth column 
presents the rms density fluctuation in an 8/i -1 Mpc top-hat 
sphere at z — 0; in the fifth column we present the power 
index n; in the sixth column we present the value of the 
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Table 1. The cosmological models chosen from Cen (1998a, b). The model numbers correspond to those presented in Cen's table 1. 
We consider a baryonic component = 0.05 and for the dark component we have = Sl m — 0.05. Models 23 and 26 include the 
cosmological constant (A). 



Model 




n A 


(T8 


n 


Ho 

(kms -1 Mpc -1 ) 


Comment 


14 


0.35 





0.80 


1.0 


70 


Open CDM 


20 


0.60 





0.73 


1.0 


60 


Open CDM 


23 


0.40 


0.60 


0.79 


0.95 


65 


A CDM 


26 


0.20 


0.80 


1.5 


0.95 


100 


A CDM 



Hubble constant; and finally the seventh column indicates 
the type of power spectrum used. 

For the normalization of the power spectra studied here 
we followed the usual methodology found in the literature, 
where the value of erg permits one to find the constant A of 
the power spectra. Our normalized spectra follow the same 
shape seen in fig. 2 of Cen (1998a). 

Thus, for a given perturbation with mass Mc, and for 
the models described in Table 1, we can obtain the present 
radius and the present density contrast. Then, to determine 
the initial conditions at the beginning of the recombination 
era, we used the general formula given by Opher, Pires & de 
Araujo (1997) for the amplification of the density perturba- 
tions (their study incorporates both growing and decaying 
modes) . 

We start the calculations at the beginning of the re- 
combination era, at redshift z lac ~ 1000, where the ioniza- 
tion degree begins to be significantly lower than unity. For 
Zi-cc > 1000 the matter is completely ionized and density 
perturbations (positive or negative) cannot evolve due to 
the photon drag that inhibits the growth of any peculiar 
velocity. We start our calculations assuming that the under- 
dense regions evolve initially with a velocity field given by 
the Hubble flow. 

In the present study we assume that initially 

Si = tfpt = Sp* = ^ 
pt pb pd 
where p t is the total density (dark matter and baryons), pb 
is the baryon density, pd is the dark matter density, and Sp t , 
8pb and 5pd are the respective density perturbations. 

The baryonic mass contained within the underdense re- 
gion is obtained from 



M h = M c n b /n n 



(30) 



and, analogously, the collisionless dark matter within the 
underdense region is 

M d = Mclid/llm. (31) 
The initial radius of the underdense region is given by 

1/3 

(l + z rcc )-\ (32) 



3 M c 



1 



4tt po (l + <5i) 



With such an initial profile, the density at the void bound- 
ary does not fall smoothly to the density of the Universe. 
Thus, the interface (between the external medium and the 
last shell of the computational grid) is a membrane-like one. 
This membrane is maintained throughout the void time evo- 
lution, because the mass inside the void region is taken to 
be constant throughout the void evolution. 



Table 2. Principal results for model 14 described in Table 1. The 
final diameter of the baryonic component is Dbf, while Ddf is 
the final diameter of the dark matter component. The density 
contrasts for the baryonic and dark matter components are, re- 
spectively, <5bf and <5fjf- 



M C (M Q ) 


£>BF(A -1 Mpc) 




D DF (/i- 1 Mpc) 


<5df 


10 12 


3.82 


-0.755 


3.89 


-0.755 


10 13 


7.18 


-0.619 


7.28 


-0.619 


10 14 


14.2 


-0.453 


15.7 


-0.455 


Table 3. 


Principal results for model 20 described in Table 1 


See 


Table 2 for notation. 








Af c (M Q ) 


D B F(fe" 1 Mpc) 


<W 


D DF (/i- 1 Mpc) 


Sdf 


10 12 


3.07 


-0.770 


3.11 


-0.770 


10 13 


5.72 


-0.609 


5.80 


-0.609 


10 14 


11.3 


-0.360 


11.5 


-0.360 



In Tables 2-5 we present the main results for the models 
described in Table 1. We see that the spectra of perturba- 
tions used here produce different results concerning the final 
diameter of the voids (and for their density contrasts) for the 
same mass perturbation Mc- 

In particular, if we consider that voids are regions with 
density contrasts S < —0.6, then large voids with diameters 
of D > 20 - 30ft _1 Mpc [the sizes found by El- Ad et al. 
(1996) and Miiller et al. (2000)] could be produced only for 
the spectrum 26 or similar one (see Tables 1 and 5). 

In fact, it is worth mentioning that even if (2 — 3)<j 
negative density perturbations are considered for the mod- 
els 14, 20 and 23 or similar ones, large void regions can- 
not be accounted for. For example, for model 23 (see Ta- 
bles 1 and 4) a 2a negative density perturbation of 1O 14 M0 
would produce a void with baryonic (dark matter) diame- 
ter of 14.7/i~ 1 Mpc (16.7/i _1 Mpc) with <5 BF = -0.639 (5 D f = 
—0.639). If instead a 3a negative perturbation is considered, 
a void with baryonic (dark matter) diameter of 17.2/i _1 Mpc 



Table 4. Principal results for model 23 described in Table 1. See 
Table 2 for notation. 



M C (M Q ) 


■DbfCi -1 Mpc) 


<5bf 


D D F(fc _1 Mpc) 


<5df 


10 12 


3.80 


-0.800 


3.90 


-0.801 


10 13 


6.84 


-0.632 


6.93 


-0.632 


10 14 


13.4 


-0.369 


14.2 


-0.372 
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Table 5. Principal results for model 26 described in Table 1. See 
Table 2 for notation. 



M C (M ) 


Dbf^-'Mpc) 




£>op(ft _1 Mpc) 




10 12 


8.30 


-0.932 


8.40 


-0.934 


10 13 


13.4 


-0.853 


13.5 


-0.853 


10 14 


23.9 


-0.709 


25.0 


-0.709 


10 15 


48.7 


-0.622 


53.5 


-0.628 



(18.9/i _1 Mpc) with Sbf = -0.776 (<5 D f = -0.776) would 
be produced. Although large, these voids are not enough to 
explain those reported in the literature. 

It is worth noting that some authors define void re- 
gions as having a density of 0.2 times the mean galaxy 
density (Schmidt, Ryden & Melott 2001). The density con- 
trast in clusters of galaxies is around 1.5-3.0, which means 
that the density contrast in void regions would amount to 
5 — —(0.5 ± 0.2), which is less conservative than our defini- 
tion. Even considering that regions with S ~ —0.3 are voids, 
it is easily seen from our results that the largest voids in the 
Universe cannot be accounted for by models 14, 20 and 23 
or similar ones. 

A possible explanation for the largest voids could be 
related to the interaction and eventual merger of two or more 
small voids. We refer the reader to a paper by Dubinski et 
al. (1993), who investigated the evolution, interaction and 
merger of voids formed from negative density perturbations. 
Then, the merger of small voids could be a way to build up 
large voids from the spectra of density perturbations related 
to models 14, 20 and 23. 

Still considering our calculations, note that, using model 
26 with 3<r perturbations for 1O 15 M0, void regions with 
diameters of up to 62/t~ 1 Mpc (baryonic component) with 
density contrast <5bf = —0.756 could be produced. In this 
case, the dark matter void would have a diameter of Ddf = 
69/i _1 Mpc (with S DF = -0.765). This result is roughly that 
inferred to the diameter of the Bootes void. 

An interesting issue has to do with the fact that the 
post-recombination Jeans mass is ~ 1O 6 M0 at recombina- 
tion and it decreases later on, and that the mass scales stud- 
ied here are in the range of 1O 12 M0 to 1O 15 M0. Thus, one 
could argue that collapse and eventual fragmentation should 
take place in the evolving (expanding) shells. To address this 
issue properly one should study in detail how density per- 
turbations evolve in the expanding shells. 

Note that it is not enough that a perturbation be Jeans 
unstable to guarantee that it will collapse and fragment. 
Consider the following example. Let us think of the evolution 
of primordial clouds after the recombination era (that is, the 
evolution of positive density perturbations). Even the clouds 
that are already Jeans unstable at the recombination era will 
initially expand before collapsing. For a cloud that stops 
expanding, its mean density (p c ) relative to the background 
is given by (see, e.g., Peebles 1993; Coles & Lucchin 1995) 



■(t) 

which corresponds to 



5.6, 



4.6 



(33) 



(34) 



This result is a lower limit since it depends only on the 
expansion rate of the Universe (the above result is taken 
for an Einstein-de Sitter universe) and on the gravity of 
the cloud. For the cases in which other physical processes 
are relevant for the evolution of the cloud, the turnaround 
should occurr for 



> 4.6. 



(35) 



In particular, we refer the reader to the papers by Oliveira 
et al. (1998a,b) who studied the collapse of Population III 
objects and consider such issues in detail. 

Concerning the void regions, their density contrast is 
always negative and they never stop expanding. As a result 
the condition given by equation (35) is never fulfilled and so 
the evolving shells do not collapse. 

Density perturbations that could be present in the 
evolving shells, however, could in principle collapse and frag- 
ment if they had the time to stop expanding and evolve. 
Note that this issue concerning the collapse and eventual 
fragmentation in the evolving shells is by itself so interest- 
ing that it deserves to be investigated in detail. We leave, 
therefore, such an issue to be considered in another paper 
to appear elsewhere. 

An interesting result concerning the distribution of dark 
matter and baryonic matter can be seen in Tables 2-5. The 
final diameter of the dark matter is always greater than 
the final diameter of the baryonic matter, when all physical 
processes considered in the present paper are included in the 
calculations. 

To see the influence of the several physical processes on 
the evolution of a given negative density perturbation, in 
particular on the final diameters of the voids, we performed 
some calculations studying individually the effects of all the 
four cooling-heating mechanisms and also the influence of 
the photon drag (Table 6). 

Note that the dissipative baryonic processes do not de- 
fine a particular characteristic void scale; their main effect 
is to segregate the dark matter and baryonic matter. 

From Tables 4 and 6 we note that the baryonic diameter 
is a little bit lower than the dark matter diameter when all 
physical processes are included, and this result remains when 
we disregard each one of the cooling-heating mechanisms. 

Another interesting result is obtained when we disre- 
gard all the processes. In this case, the final diameter of 
the baryonic matter is almost the same as the dark matter 
component. When we disregard all the processes and the 
pressure of the Universe, the final diameters are the same 
for both components, as expected, since the two fluids in 
these circumstances are subject only to gravity. 

In Table 6 we see the influence of each physical pro- 
cess one at a time. When we disregard the photon drag, 
but maintain all the other processes and the pressure of the 
Universe, the final radius of the baryonic (dark) component 
goes from 6.84 (6.94) /i _1 Mpc to 6.93 (6.99) /i _1 Mpc. This 
occurs because the photon drag acts against the expansion 
of the negative density perturbations, inhibiting the growth 
of peculiar expansion velocities. 

The cooling mechanisms Lb., L a and Lh2 have very sim- 
ilar effects on the void evolution, and their principal influ- 
ence is related to the thermal evolution of the matter inside 
the void regions. The Compton heating-cooling when disre- 
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Table 6. The influence of the physical processes on the evolution of the final diameter and density contrast of a void formed from 
a perturbation with mass 1O 13 M0 in model 23. The cooling-heating mechanisms for the baryonic matter are: the cooling from 
recombination, L a the Lyman a cooling, Lq the Compton cooling-heating, and Ljj2 the cooling by molecular hydrogen. In the second 
and third columns we present the results without all the physical processes. The superscript +U (— U) in these columns means that we 
keep (remove) the pressure of the Universe acting on the 'void boundary'. 







Without all 
processes (+ u ) 


Without all 
processes ( _u ) 


Without 
photon drag 


Without 

Lb. 


Without 
L a 


Without 
L C 


Without 

Lh2 


D BF (h- 


»Mpc) 


6.97 


7.03 


6.94 


6.86 


6.86 


6.83 


6.86 


<5bf 




-0.647 


-0.648 


-0.643 


-0.640 


-0.640 


-0.639 


-0.640 


D DF (h- 


^Mpc) 


7.02 


7.03 


6.99 


6.93 


6.93 


6.90 


6.93 






-0.648 


-0.648 


-0.643 


-0.641 


-0.641 


-0.640 


-0.641 



garded reduces slightly the dimension of the void, because 
at high redshifts this physical process maintains the tem- 
perature of the matter inside the void close to the radiation 
temperature. 

The above results also appear when we perform the 
present analysis for other models described in Table 1. In 
particular, independent of the spectra used, the physical 
processes produce very similar results concerning the final 
diameters of the baryonic and dark matter. We note that 
even using other spectra, like those used by de Araujo & 
Opher (1993) for example, this difference in diameter due to 
the presence of physical processes is maintained. 

With all the physical processes included, there exists 
a transition region with thickness of up to ~ 2.5 ft _1 Mpc, 
which depends on the mass of the perturbation, the normal- 
ization of the spectrum and the respective density parame- 
ters for the dark and baryonic components as well as the cos- 
mological constant. These transition regions have sizes com- 
parable to the diameters of galaxies and clusters of galaxies. 

The model 26, for example, has a transition region 
of ~ 2.5/i _1 Mpc thickness for a perturbation with mass 
1O 15 M . Within this region the density contrast of the cold 
dark matter is <5df = —0.628 (its diameter is 53.5 fe _1 Mpc). 
Thus, the mass density of the dark matter is pd = 0.372pd u , 
where pd u is the ambient mass density of the dark matter 
(that is, the mass density of the dark matter present in the 
Universe) . 

As the boundary of the baryonic matter is at 
48.7 h~ Mpc, within the transition region the density con- 
trast of the baryonic matter is 5b = 0, and so pb = pbu- The 
relation between the baryonic and dark components present 
within this transition region is given by 

P± = 2.70 x (^A =2.70 x = 0.90 (36) 

pd \Pdu J \\ldJ 

Thus, in this region the amount of baryonic matter is 
approximately equal to the amount of non-baryonic dark 
matter, although this universe model has three times more 
dark matter (fid = 0.15) than baryonic matter (fib = 0.05). 

Consequently, galaxies formed near voids and within 
the transition region could, during their formation process, 
have a lower content of cold dark matter than those galaxies 
formed in regions far from voids. This result leads us to 
conclude that, if one uses an analysis of the content of dark 
to baryonic matter for the galaxies in different regions, one 
could infer very different results, reflecting therefore only 
local values. 

Let us consider another example. The transition zone 



of model 26, for a perturbation with mass 10 Mq, gives 
Pb/pd ~ 7(fib/fid) ~ 2.3, although in such a universe model 
the amount of dark matter is three times greater than the 
amount of baryonic matter. 

It is worth stressing that galaxies formed out of the 
transition zone and that escape to there could have a dif- 
ferent relation between the baryonic and dark components 
than that found in the transition zone. 

It is important to have in mind that the process of struc- 
ture formation of the Universe is a very complicated one. 
For the canonical model of cosmology, namely primordial 
Gaussian random perturbation field, the formation of voids 
and clusters of galaxies is part of one process of clustering. 
The formation of regions with positive density contrasts, like 
clusters of galaxies, contributes to the growth of underdense 
regions and vice versa. 

The present study only considers a small part of the 
very complicated structure formation process. We study the 
evolution of negative density perturbations taking into ac- 
count all relevant physical processes acting on them. How- 
ever, the assessment of the role of these very physical pro- 
cesses not considered in the literature, besides including the 
presence of non-baryonic dark matter and cosmological con- 
stant, are the main contributions of the present study. 



4 CONCLUSIONS 

In the present investigation, we study the evolution of neg- 
ative density perturbations taking into account a series of 
physical processes that are present during and after the re- 
combination era. In order to perform such a study, a spher- 
ical Lagrangian hydrodynamical code has been written to 
follow the evolution of these density perturbations. 

We analysed a set of C05£^normalized spectra ob- 
tained from a recent study performed by Cen (1998a,b). The 
results presented here show that negative density perturba- 
tions can directly create large voids, with diameters of up to 
~ 50 /i Mpc, only for some universe models dominated by 
the cosmological constant (Qa ~ 0.8). 

Empty regions as large as the Bootes void could have 
been formed only from model 26 with 3<r perturbations (for 
a mass perturbation of 1O 15 M0). The diameter of the bary- 
onic component, Dbf = 62/i _1 Mpc, and the depth obtained 
from our simulations (Sbf = —0.76) are similar to those 
determined by Dey et al.(1990) for that large void. 

Other CDM spectra studied here cannot account for 
the large voids found by El- Ad et al.(1996) and Muller et 
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al. (2000) in their studies. In particular, models 14, 20 and 
23 could not account for the large voids even considering 
(2 — 3)<t negative density perturbations. The formation of 
the largest voids could come, for these spectra, only from 
mergers of two or more small voids as addressed, for exam- 
ple, by Dubinski et al. (1993). 

Our models show that the dark matter void is greater 
than the baryonic matter void, this result being produced 
by the physical processes acting on the baryons. We find a 
difference between the radii of ~ 1 — 10 per cent, depending 
on the power spectra and masses of the perturbations, when 
all physical processes are included. This effect corresponds 
to the creation of a transition zone defined by the radii of 
the baryonic and dark components. We note that disregard- 
ing each one of the physical processes, the thickness of this 
transition zone decreases and it goes to zero when all pro- 
cesses and the pressure of the Universe are removed. The 
analysis of the physical processes lead us to conclude that 
the principal mechanism that acts on the evolution of voids 
is photon drag. When we disregard this mechanism the final 
diameter of the baryonic matter increases and the thickness 
of the transition zone becomes almost half of its initial value. 

Although the thickness of the transition zone is small 
(in the range of ~ 0.03 — 2.5/i -1 Mpc) when compared to 
the radii of the baryonic and dark components, interesting 
effects on the formation of galaxies could occur. Putative 
objects formed near voids and within the transition region 
would have a dark matter to baryonic matter ratio such that 
Pd/pb 7^ fld/f2b- This means that the amount of dark matter 
is different in such regions as compared to that present in 
the Universe as a whole. 

Then, if one uses these galaxies to determine by dynam- 
ical effects or other techniques the quantity of dark matter 
present in the Universe, the result obtained would be only lo- 
cal and not representative of the Universe as a whole. Note, 
however, that within voids the relation between baryonic 
matter and dark matter maintains the same value as that of 
the Universe. 

For different spectra one would obtain different results 
as compared to those present here, but void regions with 
diameters of tens of megaparsecs and transition zones of 
non-negligible sizes would occur as well. 

Owing to the fact that the baryonic matter undergoes 
a series of physical processes, such as those considered here, 
that the non-baryonic dark matter does not undergo, mo- 
tions of matter on large scale due to the presence of a grav- 
itational potential could induce a segregation between the 
two components, producing regions with different values of 
Pb / pd that are different from the Sib / f^d relation of the Uni- 
verse. 

It is important to note that the voids studied here are 
formed only by the evolution of negative density perturba- 
tions. The voids formed in this scenario are, therefore, cold. 
We do not consider a possible energy injection, or excess 
of ionization (see, e.g., Srianand 1997), by quasars or other 
primordial objects near the voids. In this way, the temper- 
atures of the matter in our model are never greater than 
the temperature of the radiation. The inclusion of an en- 
ergy source, such as quasars, is an interesting possibility for 
a future extension of the scenario presented here. 

Certainly the process of formation of the structures of 
the Universe is a very complicated one. As discussed above, 



the formation of voids and clusters of galaxies is part of the 
same process of clustering. The formation of regions with 
positive density contrasts, like galaxies and clusters of galax- 
ies, contributes to the growth of underdense regions. On the 
other hand, the underdense regions push the matter around 
them during their evolution and growth, contributing to the 
clustering processes. 

In our study we only consider the evolution of nega- 
tive density perturbations, which is a small part of the very 
complicated structure formation process. However, the as- 
sessment of the role of a series of physical processes not 
considered in the literature, besides including the presence 
of non-baryonic dark matter and cosmological constant, are 
the main contributions of the present study. 

Also, because our study shows that the physical pro- 
cesses considered here are relevant to the evolution of nega- 
tive density perturbations, this leads us to conclude that in 
the structure formation of the Universe as a whole these very 
processes could be relevant as well. Therefore, it would be 
of interest to take into account, in more realistic modelling 
of the structure formation, the various processes considered 
here. 
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